Amplification of compressional MHD waves in systems with forced entropy oscillations 
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The propagation of compressional MHD waves is studied for an externally driven system. It 
is assumed that the combined action of the external sources and sinks of the entropy results in 
the harmonic oscillation of the entropy (and temperature) in the system. It is found that with 
the appropriate resonant conditions fast and slow waves get amplified due to the phenomenon of 
parametric resonance. Besides, it is shown that the considered waves are mutually coupled as a 
consequence of the nonequilibrium state of the background medium. The coupling is strongest 
when the plasma /3 ~ 1. The proposed formalism is sufflciently general and can be applied for many 
dynamical systems, both under terrestrial and astrophysical conditions. 
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I. INTRODUCTION 

Entropy plays a central role in the thermodynamic 
properties of fluids and plasmas that consist of an enor- 
mous number of particles. In general terms and as a 
function of energy, volume, etc., the entropy is a mea- 
sure of disorder as it counts the multitude of microscopic 
realizations for some given manifest condition. That dis- 
order is realized in the chaotic motions of the system 
entities and in the dissipative character of transport pro- 
cesses. Under normal circumstances, vifhen the system is 
left to itself, it evolves to the state of maximum entropy 
which is compatible with the imposed conditions. That 
relaxation to a so-called equilibrium is one realization of 
the direction of time and it could be applied to any closed 
system or to the universe as a whole for that matter. 

Open systems, however, can be found and for a long 
time in a state of low entropy. In that case, these systems 
reach a steady state regime in which they are constantly 
driven out of equilibrium. An important topic for investi- 
gation is the study of various dissipative effects and their 
relation with fluctuation and transport phenomena. It is 
of great interest to understand if and how the low entropy 
is realized in some type of macroscopic pattern or struc- 
ture, i.e. a collective and ordered motion which stands in 
great contrast with the chaotic motion at a microscopic 
level. Obviously, macroscopic order is also possible at 
equilibrium, e.g. at low temperature or at high pressure, 
but in the present paper we aim to describe a nonequi- 
librium effect in which a nontrivial macroscopic effect is 
generated in time. More specificly, our study involves the 
ampliflcation of wave motions in a driven system. 

In order to organize the chaotic motions of the particles 
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into collective motions exhibiting long range or long term 
non-trivial correlations, the effort of some external driver 
is required. This external source supports the system by 
delivering work or exchanging heat. Secondly, the system 
must react to the driver in such a way that the interac- 
tion between the different structures is sufficiently rich 
to select a nontrivial mode of operation. The specifics 
needed of that nonlinear dynamics remain more vague as 
today there is no unifying understanding of morphogene- 
sis. In the present paper, the mechanism at work is that 
of parametric resonance, as we will explain below. 

We consider compressional magnetohydrodynamic 
(MHD) waves (structures) propagating in a magnetized 
plasma (the system consisting of chaotically moving ions 
and electrons), which is driven by external heat reser- 
voir(s). Mechanisms of excitation, conversion, ampli- 
fication etc. of different wave modes are of special im- 
portance for many terrestrial and astrophysical systems. 
The waves represent effective transmitters of the energy 
and of the information present in the background sys- 
tem. There are many studies concerning the wave prop- 
agation in driven systems. One of the many examples 
is shear flow. In a shear flow, an inhomogeneous veloc- 
ity proflle is formed as a consequence of the action of 
the external force. The flow is steady and the system is 
in a stationary nonequilibrium state. In particular, the 
work done by the external force does not increase the 
system's energy. On the contrary, the delivered energy is 
dissipated and, more interestingly, different wave modes 
exchange energy with each other and are coupled by the 
flow , see for instance [1, Q and references therein, and 
d, [23, [2I] for more applications in the context of the 
nonmodal and pseudospectral analysis of such systems; 
the role of shear flows in the formation of the observed 
p-mode power spectra is discussed in [20] and its role in 
the heating of the magnetically governed solar corona is 
found in [22 |. 

One of the essential mechanical properties of fluids 
and plasmas is the possibility of parametric resonance. 
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which might occur in a system sustaining different kind 
of periodic in time (or space) processes. Recently, the 
abihty of different types of MHD waves to be involved 
in so-called swin g (o r parametric) interactions, has been 
studied e.g. in l33l [341 ■ including spatially inhomoge- 
neous backgrounds [2ll |. However, in all these cases the 
problem was concerned with non-driven systems, with 
waves representing the perturbations above some equi- 
librium configurations of the plasma. When an external 
periodic force on the compressional fast and slow magne- 
tosonic waves was considered [s^ , nevertheless the exter- 
nal force applied to the system was treated as adiabatic, 
i.e. , it did not make any contribution into the variation 
of entropy (or temperature). Therefore, even with this 
action the characteristic sound speed remained constant 
and the parametric action of this force on the waves had 
been achieved by the periodic oscillation of the modal 
wavelength only. 

In the present paper we want to understand and to 
clarify a similar process when the external driver is a 
cause of entropy (or temperature) variation. Some ef- 
forts have already been taken to extend the model in 
this direction before. For instance, studies have been 
conducted on the parametric amplification of acoustic 
waves [H [l^, also including systems that are exposed 
to periodic variations of external electric (piezoelectric 
medium) [l^ or magnetic (magnetostrictive solids) [2§| 
fields. These variations indeed give birth to oscillations of 
the phase speeds of the waves leading to the parametric 
amplification of them under certain resonant conditions. 
In spite of these investigations, a more general theory of 
parametric amplification of fast and slow magneto-sonic 
waves in systems with periodic entropy production has 
not been developed. We address this issue in a systematic 
manner to develop a unified formalism for the qualitative 
and quantitative treatment of the problem. 

The paper is organized in the following way: in the 
next section HI] we develop a basic MHD model and 
specify the properties of the background nonequilibrium 
state; Further, in section |TTT] we study the propagation of 
fast and slow magnetosonic waves in the system (also we 
give a discussions about the hydrodynamic limit); And 
in section HVl we give some conclusive remarks. 



II. BASIC MODEL 

The aim of the present paper is to study properties of 
the compressional fast and slow magneto-sonic waves in 
a magnetized single component plasma, which is driven 
by one or more external reservoirs that bring it out of 
thermodynamic equilibrium. We can write the full set of 
MHD equations to describe the dynamics of such plasmas 
in the following way: 
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where all the quantities have their usual meaning and 
where D/Dt = d/dt + (V ■ V) denotes the convective 
derivative. The symbols v and rj denote the coefficients of 
viscosity and magnetic diffusion, respectively. However, 
we write these dissipative terms in the equations just 
for the sake of generality; as will be shown below, the 
effects of viscosity and resistivity can be omitted in this 
consideration. Such a more general set-up is necessary to 
indicate that, in general, the mentioned effects can also 
give rise to a contribution into the source term, C (i.e., in 
the production of entropy), standing on the right-hand 
side (RHS) of Eq. ^. 



A. The background state 

Before we turn to the study of the wave properties in 
the considered nonequilibrium system, it is convenient 
to specify the background state. In fact, this is desir- 
able because usually, in the literature, the wave pro- 
cesses are considered to be adiabatic, i.e. the action of 
non-adiabatic processes, leading to the generation and/or 
transport of the entropy within the system and its sur- 
roundings, is usually omitted. This approach implies a 
definite separation between two characteristic timescales, 
viz. (i) that of the non-adiabatic variation of the entropy 
and (ii) the timescale of collective (macroscopic, mechani- 
cal) phenomena, for instance the wave and the oscillatory 
motions. From the structure of the basic set of equations 
(IT|)-(|3]), it is evident that we consider also non-adiabatic 
processes, which can bring the system away from the 
thermodynamical equilibrium. In general, Eq. ([¥]) rep- 
resents an entropy balance equation of the form: 
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where T and S denote the temperature and entropy per 
unit volume, respectively. The source term in this equa- 
tion represents the sum of all the sources and sinks of the 
local entropy. As is well known from the characteristic 
analysis, when one considers only the adiabatic motions 
of a plasma or a fluid (£ ~ 0, or in other words, when the 
background state is in exact equilibrium) then any small 
perturbations of the entropy are aperiodic (often referred 
to as entropy waves/modes) and they are decoupled from 
the other fundamental eigenmodes of the system (which 
are of oscillatory nature). Hence, in that case, the local 
entropy remains unchanged as these other eigenmodes 
propagate. 
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When the system is open, however, the local processes 
can be considered as part of a larger system which itself is 
closed but still enables local variations (even decreasing) 
of the entropy. When the characteristic timescales of the 
considered processes are comparable with the timescale 
of the entropy variation, the source term in Eq. ((5]) can 
be significant and, thus, should not be neglected in that 
case. In this situation, the entropy variations can be 
coupled with the other modes of collective motion which 
are sustained by the system. If the temporal variation of 
the entropy source is harmonic, then one can refer to the 
entropy deviations from its equilibrium value as 'forced 
entropy modes' (or oscillations). 

In general, one can consider a static equilibrium char- 
acterized by the entropy Sqq , with an homogeneous pres- 
sure Poo = const and density poo = const. In this 
equilibrium, the source term vanishes, i.e. £00 = 0. 
Further, we introduce a small time-dependent deviation 
from the equilibrium, so that 5*0 = Sqq + Soi{t). This 
deviation is due to the presence of external reservoirs 
which exchange entropy with the system leading to a fi- 
nite time-dependent perturbation modelled by the term 
£01 (t). The physical quantities po = poo + Poi(t) and 
Po = Poo + Poi(i) then satisfy the following equation: 
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This equation, in combination with the appropriate equa- 
tions from the basic set (H])-©, implies that the entropy 
variations in general are not decoupled from the motions 
related to the fundamental modes in the system. More- 
over, under certain conditions (which we will specify be- 
low) these equations can describe forced oscillations. [38] 
Now let us consider the source term Coi{t) in more 
detail. There are different (transport) processes which 
may contribute to this term, which formally could be 
attributed to two classes, viz. (i) the processes of the en- 
tropy production (such as viscous and resistive dissipa- 
tion, thermonuclear reactions, chemicals reactions etc.), 
and (ii) the transport of entropy (via thermal conduc- 
tion, radiation, diffusion, etc.). Following the represen- 
tation given by Groot and Mazur [3], this term can be 
represented in its most general form as: 



£o(i) =To(-(V- J,o)+^To) 



(7) 



Here, J^o denotes the entropy current, and cto denotes 
the entropy production. We should note that 3s — 
Js.tot — pSY, i.e. 3s represents only the microscopic cur- 
rent of the entropy, while the macroscopic transport pS\ 
is excluded [7]. 

Further, we represent both these quantities as the sum 
of a constant and a variable part, i.e. ctq = ctoo + croi{t), 
and 3so — Jsoo + Jsoi(l)- These are the general rep- 
resentations of the entropy production and the entropy 
current. In this paper, we are interested to study systems 
in which the constant parts of these two effects compen- 
sate each other, supporting the development of thermo- 
dynamic equilibrium when entropy variations in time are 
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and the remaining part of the source term, £01 j depends 
on time in a harmonic way: 

£01 =To(-(V- J,oi)+^Toi) --asinm, (9) 

where is the frequency of the entropy variation in time 
and a denotes a constant amplitude. This amplitude is 
rather small as we consider only small oscillations around 
thermodynamic equilibrium. Naturally, a might depend 
also on the spatial coordinates leading to a spatial varia- 
tion of the thermodynamic quantities as well. However, 
we here assume that this variation is rather weak. In 
other words, the characteristic length of this variation is 
assumed to be larger than the system size and the action 
of the external reservoirs thus causes an oscillation of the 
entropy everywhere within the system at once. Conse- 
quently, in the equations governing the dynamics of the 
background state, we replace the convective derivatives 
by partial derivatives. 

Now, Equation ^ reveals the fact that the considered 
variation of the entropy can result in a variation of the 
other physical quantities, for instance the plasma pres- 
sure, density, temperature, velocity, etc. The purpose 
of the present work is to discover the basic nature of the 
resonant processes which may take place in the described 
dynamical system. And, therefore, and for the sake of 
further simplification, we assume that the entropy oscil- 
lation leads only to an oscillation of the pressure and the 
temperature, while the plasma density remains constant 
(i.e. poiit) = 0). This is possible in any configuration, 
which admits only a very limited ability for the system 
volume to vary significantly or when the system is not 
able to alter its volume at all globally in the background 
(isochoric process), while small fluctuations of the den- 
sity are still possible. With this assumption, Eq. ^ takes 
the form: 



dt 



= —ai sin{nt), 



(10) 



where, ai = (7 — l)a. This equation can be integrated 
to obtain pq = pqq +p» cos{ilt) (here, = ai/fl), which 
immediately implies that the characteristic sound speed 
of the system is a harmonic function of time as well 
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where Csq and S are constant parameters. 

One might be interested in an example of such a con- 
figuration where this kind of situation can be realized. 
Consider, for instance, a static (Uq = 0) system with 
thermal conduction driven by the heat flux through the 
boundaries. For such a system one can then write down 
the following Fourier's equation: 
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where C„ denotes the specific heat at constant volume 
(density), A is the coefficient of the thermal conduction, 
and A denotes the Laplacian This equation is one 
realization of Eq. In the appendix, we also show 

the relation between the coefficient of thermal conduction 
and the local production and transport of entropy. So, if 
the boundary condition is a harmonic function of time, 
then the temperature profile within the system would be 
also oscillating representing an example of the realization 
of a set-up given by Eq. pUj) . 

Consider the case when the temperature at the bound- 
ary a; = is a periodic fimction of time: 



Tq — Too + Toie 



0. 



(13) 



The distribution of the temperature within the system 
then has the same time dependence, i.e. exp —iujt. The 
real temperature distribution will then be [cf. 
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Therefore, it is seen that the temperature oscillation 
propagates in the form of damped thermal waves. Here, 
one can avoid the spatial dependence in a similar way 
as it was done by Lorenz to stud y th e properties of the 
Rayleigh-Benard convection cells [1J|. However, this can 
be done only when the wavelength of the fundamental 
harmonic of the temperature waves given by Eq. (|14p is 
larger than the system size (i.e. when either the coeffi- 
cient A is large or the density (or oscillation frequency) 
is small enough). 

There is also another way to achieve the variability 
of the background temperature: if we have an external 
source of heat Q (for example, in the laboratory situa- 
tion, this source can be the heating by electric currents 
or lasers), then Eq. p2|) can be rewritten as 



dT 

pC,— ^XAT + Q. 
at 



(15) 



Now, if we take Q as periodic function of time and homo- 
geneous in space, then the temperature (pressure) may 
have a periodic time dependence as well, thus leading to 
Eqs. (Uni) and dU). 



III. FAST AND SLOW MAGNETO-SONIC 
WAVES 

In this section, we derive the dynamical equations gov- 
erning the propagation of the linear eigenwave modes of 
the system. For this purpose we linearize the full set 
of MHD equations over the background nonequilibrium 
state outlined in the previous section. In principle, the 
process we are addressing here represents a weakly non- 
linear action of the entropy oscillations, forced by the 
combined effect of all sources and sinks of entropy upon 
the eigenmodes of the system. At this stage, we treat 



the problem in two spatial dimensions only to avoid un- 
necessary mathematical complications. The plasma is 
embedded into a uniform magnetic field directed in the 
z-direction. The x-axis is assumed to cover the direction 
across the magnetic field. We know that, in this partic- 
ular case, the system sustains only two kinds of funda- 
mental modes, viz. fast and slow magneto-sonic waves. 
The linearized set of equations then reads: 
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Po dt 
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+ vAux, (17) 
(18) 
(19) 
(20) 
(21) 



where the perturbed vector fields of velocity and mag- 
netic field are written by lowercase symbols, while the 
scalar fields (pressure and density) are denoted by sym- 
bols with primes. Further, we assume that the wave- 
length of the waves that are considered here are larger 
than the characteristic dissipation scales in the system. 
Therefore, we can neglect the terms corresponding to 
the viscosity vAux , vAu, ~ and the magnetic diffu- 
sion riAbx,riAbz ~ 0. Consequently, we assume that the 
perturbation of the source term in Eq. (|2ip vanishes, i.e. 
C w 0, as it represents a contribution in the entropy vari- 
ation from the dissipation of the perturbations. In other 
words, we consider the situation when the energy of the 
external reservoirs driving the system is much larger than 
the energy of the waves and a back reaction of the lat- 
ter on the system is negligibly small. Yet, the effect of 
the variable background is represented by the variable 
pressure or sound speed. 

The governing set of equations (|16 p - (PT1) is homoge- 
neous w. r. t. the spatial coordinates. Hence, we can ap- 
ply a Fourier analysis by representing all physical quan- 
tities as: 



/ = 



f{kx, kz,t) expi{kxX + kzz)dkxdkz, (22) 



and with the above mentioned assumptions we obtain a 
set of two second order ODEs: 



+ CiiUx + Ci2Uz = 0, 



(23) 
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+ Cl2Ux + C22Uz = 0, 



(24) 



where cn — C'^{t)kl 
Cg{t)k'^ (note that here, 
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kr,. 



k1 and Va 

Bo/y/ATTpo corresponds to the Alfven speed). These 
equations govern the evolution of the osciUatory motions. 
Before we make a rigorous analysis of this system we con- 
sider the hydrodynamic case, i.e. when Bq = 0. 



A. The hydrodynamic case 

When the magnetic field is absent, the system 
can be rewritten in the form of a single wave equation. 
Introducing the new variable: 



kxUx ~t~ kzUxi 



we can write: 



^ac + c^'ofc'(i + Scos{nt))x,, = 0, 



(25) 



(26) 



appropriate eigenvectors can be constructed which cor- 
respond to eigenvalues calculated from the characteris- 
tic equation and which are perpendicular to each other. 
Using this procedure, one can then show that the funda- 
mental modes are separated (or decoupled). 

However, when the background system is driven (i.e. 
when S ^ 0), the coefficients in the equations are time 
dependent. Therefore, the analysis of these equations 
becomes rather complicated in this case. We are going 
to construct the 'eigenfunctions' by employing a similar 
procedure in combination with numerical techniques, in 
order to draw a more complete picture of the dynam- 
ical processes we are studying here. For this purpose 
it is convenient to impose an appropriate normalization 
of the quantities. We introduce the following scaling of 
the variables: = u^/Va, Uz = Uz/Va, t — VAkt, 
Cii = cn/VAfc = + 1, Ci2 = ci2lVAk = eK^K,, 

C22 = C22/VAk = S.'^K'^, here, Ki = h/k {i = x,z), 
^2 = C^/V^ = ^2(1 + 5 cos(f]T)), and n = n/VAk. Ob- 
serve that, if we apply an orthogonal transformation on 
the velocity vector space {Ux, Uz) of the form 0, [ll|: 



where the number of dots above the variable indicates 
the order of time-derivative. This reduction to a sin- 
gle wave equation represents the fact that there are only 
compressional acoustic waves in the hydrodynamic limit. 
It is evident that Eq. ((26)) is a Mathieu-type equation 
and, thus, an exact analytical solution is available. The 
solution has a resonant nature when: 



Xf = cos 9Ux — sin 6Uz, 
Xs = sin 9Ux + cos 9Uz, 



(30) 



where 9 denotes the Euler angle in the velocity space. 
This transformation is time dependent since the coeffi- 
cients in the initial matrix: 



n^2Cso\k\, 



(27) 



and, more precisely, when the basic frequency of the 
acoustic wave lies within the following interval 



Cso\k\ - — 



< 



(28) 



also depend on time: 



tan 6* = 



Cii C12 
C12 Cii 



C12 



(31) 
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12 



fij ~ Cll C^22 ~ f^i 



where /i = SCg^k . With this notation, the resonant 
solution of Eq. (pS)) takes the following form: 



cos I — T 



smi-r 



Xac(i)=^ac(t-0)exp(^^T 

(29) 

As a conclusion, when the entropy (temperature) oscil- 
lates with a frequency that is twice the basic frequency 
of the acoustic wave, then those waves are resonantly 
amplified in the system^ 



B. Separation of MHD waves 

Let us now turn to the general analysis of the system 
([23|l -(|24 p . This system consists of two wave equations 
which describe the presence of two different wave modes. 
The velocity fields corresponding to these two kinds of os- 
cillatory motion of the system are clearly coupled. Nev- 
ertheless, if one considers a perturbation around a static 
equilibrium state (i.e. with J = 0), then it is known that 



Cii — rif 

C12 



~ ^^22 

C12 



(32) 



With this notation, the initial set of equations 
can be transformed into: 



X. 



nl-9'^] Xf = 9Xs + 29Xs, 



(33) 



Xs + {hi - 9^^ Xs - -9Xf - 29Xf, (34) 



where rip and Og are the eigenvalues of the matrix ([5T|) 
representing, in fact, the characteristic frequencies of the 
fast (F) and slow (S) magneto-sonic waves, respectively: 



1 



^^F,s = ^ [Cii + C22 ± V (^11 - ^11)' + 4C?2 j . (35) 

It is well known that in a magnetized plasma there ex- 
ist three different regimes, which are determined by the 



6 




FIG. 1: Xp (solid lines) and Xf (dotted lines) plotted against time r. 



relative importance of the thermal and magnetic pres- 
sures. This relation usually is parameterized by the ratio 
of the thermal and magnetic pressures, the plasma beta 
defined as /3 = Snpo/B^ = 7C^/2. The proportionality 
factor here is of the order of unity. The three regimes 
are referred to as » 1, <C 1, and ~ 1, respec- 
tively. Accordingly, the basic properties of the different 
MHD waves are also well studied in these three regimes. 
Below, we will also consider the problem case by case. 
Typically below we use instead of plasma (3. 



C. The case > 1 

Firstly, let us consider the limit of large plasma beta. 
In this case, the magnetization of the plasma is rather 
weak and the thermal effects dominate in the dynamics of 
the medium. In this regime, the fast and slow magneto- 
sonic waves have drastically different properties. In this 
situation, one obtains from the expressions (|35p that 

> (36) 



while, on the other hand, the low frequency branch of the 
spectrum satisfies: 

^ 1, (37) 

showing that the characteristic frequencies of the slow 
waves lie close to the Alfven frequency. In general, the 
fast and slow magneto-sonic modes propagating in an ar- 
bitrary, oblique direction with respect to an applied mag- 
netic field are neither purely compressible (longitudinal) 
nor incompressible (transversal) . The latter property re- 
lates to the tension of magnetic field lines, like for Alfven 
waves. But, both of the wave modes consist of a mix- 
ture of these two features. However, when ^ V^, the 
properties of the fast magneto-sonic waves are very sim- 
ilar to those of sound waves and these waves represent 
sound waves that are slightly modified by the presence of 
the magnetic field and the transversal component is weak. 
Therefore, these waves substantially feel the variation of 
the sound speed, what directly results in a variable phase 
speed of the fast modes. 

We now analyze the set of equations (P5)) - (|34p term by 
term. For this purpose one should use (5 as a small param- 
eter in the problem, as we consider here a nonequilibrium 
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FIG. 2: As in Fig. [T] for ^ < 1. 



background that is very close to the equilibrium (5 <C 1. 
We are interested in terms of the order of magnitude up 
to linear in i5. All higher order terms yield vanishingly 
small contributions to the considered dynamical process. 
It can be easily shown that includes both zeroth- and 
first-order terms in the expansion in terms of 6, while the 
terms 9 and 6 do not contain zeroth-order terms and the 
expansions of both these terms start with the first-order 
term. Consequently, the term 9^ can be neglected in both 
equations as it consists of terms definitely of second and 
higher orders in S. Further, we assume that within the 
resonance, the fast magneto-sonic wave would be able to 
grow in amplitude exponentially, while the terms con- 
trolling the coupling with the slow waves (on the RHS) 
are rather small as the characteristic frequencies of these 
eigenmodes lie far apart from each other according to 
Eqs. ([55]) - ([57]) . Consequently, we neglect the couphng 
terms, which at any rate will be dominated by the terms 
on the LHS of Eq. ((33)) . and we are left with two approx- 
imate 'decoupled 'equations. With these simplifications, 
Eq. ([55]) reduces to: 



(38) 



where np — np{S,fl) <^ 1 and ip-p — ip-p(S,Q) are the 
amplitude and the initial phase of the fast wave frequency 
oscillations in the system, respectively, and ilpo is the 
equilibrium frequency of the fast magneto-sonic waves. 
This Mathieu-type equation has a resonant solution for 
the appropriate harmonic of the fast magneto-sonic wave 
with the frequency (corresponding to a given value of the 
wavelength) : 



varying within the interval 



FO 



(39) 



(40) 



and reading as: 

Xpit) =Xf{t = 0)exp 



Xf + ClpQ (1 -I- np cos(riT -I- ipp)) Xf = 0, 



n 

cos I — r 
2 




(41) 
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The approximative analysis, which is given here, 
should be verified by means of experimental or numer- 
ical methods. In the present study, we first concentrate 
on a rigorous numerical analysis of the initial equations 
(l^5)l - (|24p . and then we reconstruct from the results of the 
numerical simulations, all the relevant quantities we use 
in the analysis. 

Here, we consider two different sets for the initial con- 
ditions. First we chose the initial conditions Uxo, UxO: 

Uzo and Uzo so that, initially, only the amplitude of the 
fast magneto-sonic wave (solid line in panel Ai) is finite 
and there is no slow magneto-sonic wave in the system 
(see the dotted line in panel Ai): 

Xs(T = 0)~e, i:s(r = 0)~£, (42) 

where e « 10~^^ is the order of magnitude of the nu- 
merical error related to the machine precision. Panel Ai 
represents the case of the stationary state in the system 
((5 = 0). On the other hand, in panels Bi and Ci of the 
same figure, the curves corresponding to the parameter 
values S = 0.05 and ClpQ = n/2 « 10 are plotted. It 
is obvious that the amplitude of the fast oscillation un- 
dergoes a rapid exponential growth in this case (panel 



Bi, Fig [T]), while the slow magneto-sonic wave is out 
of resonance and, consequently, its amplitude remains 
small (panel Ci). The excitation of this small amplitude 
slow magneto-sonic disturbance is due to the weak cou- 
pling of the fast and slow waves by the external driver. 
This coupling implies that, even when there is initially 
no slow magneto-sonic wave, nevertheless, a small ampli- 
tude slow magneto-sonic wave perturbation gets excited. 
These small amplitude vibrations can be referred to as 
'slow' magneto-sonic noise. The presence of this noise re- 
lates to two simultaneous factors, viz. the propagation of 
the fast magneto-sonic wave and the external driving of 
the background system. The exclusion of either of these 
two factors would lead to the disappearance of the slow 
magneto-sonic perturbation, as is shown in panel Ai. 

On the other hand, we have carried out the simula- 
tions for the case when frequency of slow waves satisfies 
the resonant condition Cl^o — fa 1. The sound speed 
yields a negligible contribution into the frequency of the 
slow magneto-sonic waves, as from expression ([37l) it is 
evident that the latter is close to the Alfven frequency 
(which in our case is constant), and whatever variations 
of the sound speed (related to the pressure variation) oc- 
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cur in the system, the slow magneto-sonic wave frequency 
remains almost unchanged. The corresponding values of 
the increment Ss is so small that there is not a significant 
growth of the slow mode amplitude during any reason- 
able time span. This means that this variation is not felt 
by slow magneto-sonic wave modes. In the slow magneto- 
sonic waves propagating in a high-beta plasma, the in- 
compressible (transversal) component dominates and the 
properties of these waves are very similar to the purely 
incompressible Alfven waves. Therefore, one can expect 
that the slow magneto-sonic waves (and also the Alfven 
waves) do not feel the variations of the sound speed when 
the Alfven speed remains constant. This is why, in the 
case when we initially have only the slow magneto-sonic 
wave {A2 of Fig. [T]), these waves stay out of the para- 
metric interaction (panel C2), even when an appropriate 
resonant conditions is formally satisfied. In panel B2 of 
Fig. [U we plot X-p showing that the fast magneto-sonic 
wave naturally is out of resonance due to the fact that its 
frequency is out of the resonant area. Again, the pres- 
ence of the slow magneto-sonic wave of finite amplitude 
in the externally driven system leads to the excitation of 
small amplitude disturbances (panel -82), which we can 



now refer to as 'fast' magneto-sonic noise. 



To involve the almost incompressible slow magneto- 
sonic waves into the parametric interaction, one should 
not consider the case of constant density (for the para- 
metric amplification of Alfven waves with periodical vari- 
ation of medium density see [s^l). Instead of Eq. (|10p . 
the full equation Eq. ([6]) should then be employed. In 
that case, the action of the external periodic forces can 
also be taken into account. However, along with the 
effect of parametric resonance, we expect that strong 
couplings between the conipressional fast and slow and 
incompressible Alfven waves can arise. Therefore, the 
problem would then not be reducible to two spatial di- 
mensions any more, since the Alfven speed also would be 
variable and the Alfven waves would participate in the 
process as well. This issue can become the subject of 
future developments and generalizations of the present 
model. 
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D. 



The case < 1 



and it has the form: 



In the present subsection, we address the case when 
the plasma is magnetically dominated, i.e., the case 
<C FJ. In this regime, the situation is just opposite 
compared to the previous one. In particular, we know 
that now the frequency of the fast magneto-sonic wave 
satisfies 



> 1, 



(43) 



and the properties of these waves are close to those of 
the incompressible Alfven waves. Consequently, in a low- 
beta plasma the fast magneto-sonic waves almost do not 
feel the variation of the entropy. However, this is again 
true only when an external driver causes a variation of 
only the pressure while the density remains constant. In 
general, when the Eq. © holds so that the density is vari- 
able, the remark given in the previous subsection about 
the slow magneto-sonic waves and Alfven waves is now 
valid for fast magneto-sonic waves. 

Within the present framework, the fast magneto-sonic 
waves stay out of resonance. This can also be discovered 
by means of numerical simulations. In Fig. [21 we show 
the results of our simulations for the new set of parameter 
values in the same order as in Fig[TJ Panel Bi confirms 
the validity of the conclusions we have just made. 

However, one can again follow an approximative anal- 
ysis for the slow magneto-sonic waves, in a similar man- 
ner as in the previous case. Indeed, contrary to the case 
of high plasma beta, now the slow magneto-sonic modes 
are predominantly compressional (longitudinal) and their 
properties are close to those of the usual sound waves: 



^s<e 



(44) 



Therefore, following a similar logic as before, one can de- 
rive an approximate Mathieu-type equation for the slow 
waves which reads: 



nlo{l + nsCOs{nT + ^s))Xs = 0, (45) 



where ns = ns{S,fl) <C 1 and ips = (pF{5,ft), respec- 
tively, are the amplitude and the initial phase of the slow 
magneto-sonic wave frequency oscillations in the system, 
and Clso is the equilibrium frequency of the slow magneto- 
sonic waves. Again, this Mathieu-type equation has a res- 
onant solution for the appropriate harmonic of the slow 
magneto-sonic wave with the frequency (corresponding 
to a given value of the wavelength) : 



varying within the interval 



n 



FO 



< 



(46) 



(47) 



Xs{t)^Xs{T^O)exp 



so 



2n 



n 

cos I — T 



sm I — T 



(48) 



Furthermore, we verify these findings by a numerical sim- 
ulation when Ctso = ^1/2 0.1. As we show in panel 
6*2 of Fig. [21 the slow magneto-sonic waves get ampli- 
fied, while the fast magneto-sonic waves stay out of reso- 
nance (panel B2), and the coupling between these waves 
is again very weak because f2so ^ ^^fo- Observe that in 
both cases of the initial conditions (panels Ai and A2), 
we can see the excitation of the small amplitude 'slow' 
magneto-sonic (panel Ci) and 'fast' magneto-sonic (panel 
B2) noises due to the weak coupling of the waves. 



E. 



The case ^ ~ 1 



This case is of special importance in the problem con- 
sidered here. When « Vj, the thermal and magnetic 
effects are of comparable strength. Under these condi- 
tions, the properties of both waves are very similar as: 



(49) 



That is the reason why the waves can not only interact 
with the forced entropy oscillations but, in addition, they 
also may effectively couple and exchange energy among 
each other. This can be understood easily if we recall 
the presence of the terms in the RHS of Eqs. ([55)) - ([M)) . 
These equations show that an external driver can effec- 
tively couple the waves. This situation is similar to the 
case of a laminar shear flow in the system j^. As is 
well-known, the latter also can be maintained only by 
an external driving of the system. This seems to be a 
general property of driven systems: the ability to couple 
eigenmodes of the system under certain conditions. This 
claim looks apparent in the framework of the current pa- 
per. However, a complete and general proof of it can 
be given only through the general consideration of the 
fully three-dimensional case when the variability of the 
background density and the action of the external forces 
also are taken into account. Here, we should emphasize 
that the coupling of different kinds of wave modes can 
be achieved not only with the external driving, which 
harmonically oscillates in time, but in general, with any 
kind of driver giving rise to a time dependent entropy 
(temperature) in the system. 

Now, a similar approximative analysis as we have dis- 
cussed in the previous two cases, would be very 'rough'. 
In the case when the waves propagate almost along the 
magnetic field {K^ <C 1), it even becomes completely in- 
applicable. Consequently, it is convenient to treat the 
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problem in this case only numerically. Before we do this, 
it is important to emphasize again that the coupling of 
the different wave modes is possible only because of the 
presence of the external driving in the system. The sim- 
ulations show that the dynamics of both wave modes 
(when the resonant conditions are satisfied) is character- 
ized by the interplay of two different factors. On the one 
hand, both wave modes are involved in the parametric 
resonance and, on the other hand, they are effectively 
coupled to each other. 

Let us first consider the case when both fast and slow 
magneto-sonic modes propagate almost parallel to the 
magnetic field {Kr,. ^1). In this case, in [3] we ac- 
cordingly show the results of the simulations as in the 
previous subsections. It is apparent that both 'eigenfre- 
quencies' lie very close to each other. As a consequence, 
in Fig. [3] we show that, if initially we have only a fast 
magneto-sonic wave (panel Ai) and the resonant condi- 
tion is satisfied (fl — 2), the fast magneto-sonic mode 
gets amplified via the parametric action (panel Bi). Be- 
sides, the coupling of the two wave modes is very strong 
in this case and, therefore, the slow magneto-sonic wave 
is excited and also becomes involved in the amplification 
process, as this wave satisfies the resonance condition as 
well (panel Ci). Both kinds of wave modes effectively 
exchange energy and the characteristic values of their 
amplitudes are comparable to each other. The situation 
is absolutely similar, when initially we set up only a slow 
magneto-sonic wave (panel ^2)- We then observe that a 
fast magneto-sonic wave is excited (panel B2) and both 
waves are in parametric resonance with the driver as they 
exchange energy (panels B2 and C2). It is important to 
note that, in this case, the ultimate 'fate' of the process 
is more sensitive to the initial conditions as, contrary to 
the case of weak coupling, not only the frequencies are 
important but also their phases matter. 

We also investigated the case when the waves prop- 
agate in the oblique direction = 0.5. In this case, 
the frequencies of the modes become somewhat separated 
even both of them still are of same order of magnitude. 
Therefore, the coupling effect weakens substantially (as 
shown in Fig. |4]), and the parametric resonance effect 
only effectively amplifies, with the respective resonant 
conditions, in the one case only the fast magneto-sonic 
wave (panel Bi) and, in the other case, only the slow 
magneto-sonic wave (panel C2). At the same time, it 
should be noticed that again we discover the generation 
of the 'slow' (panel Ci) and 'fast' (panel B2) magneto- 
sonic noises, which are of much larger amplitude now (as 
expected) than in the limits of high or low plasma beta. 



IV. CONCLUSIONS 

In the present paper we investigated wave motions and 
wave amplifications in a non-equilibrium setting. In par- 
ticular, we have developed a model of a magnetohydrody- 
namic system that is being driven by the combined action 



of external sources and sinks of heat. More specifically, 
the external action is a harmonically oscillating function 
of time. In contrast with many other studies, the back- 
ground is that of a non-equilibrium state against which 
we have considered the dynamical properties of the com- 
pressional fast and slow magneto-sonic waves. Of spe- 
cial interest has been the case when the external driving 
of the system results only in a variation of the pressure 
(or temperature) and it does not lead to a distortion of 
the background velocity and density fields. A numerical 
study has led us to the following conclusions: 

1. In the hydrodynamic limit [Bq — 0), the spectrum 
of the possible wave modes reduces to that of the 
sound waves. In this case, we derived an exact 
equation, viz. Eq. (j26p with an obvious resonant 
solution given by Eq. 

2. In the case of a high beta plasma (^^ ^ 1), the 
temperature oscillation leads to a major contribu- 
tion in the variability of the fast magneto-sonic 
wave frequency fip, while the amplitude of the 
slow magneto-sonic frequency oscillation is very 
small. As a consequence, in this regime only the 
fast magneto-sonic waves can participate in the 
parametric resonance with the entropy oscillations 
(panel Bi in Fig. [T|) as is predicted by the approxi- 
mate solution (PT|) . The slow magneto-sonic waves 
stay out of resonance because the background den- 
sity (i.e. Alfven speed) is assumed to be constant 
(cf. panel C2 in Fig. [T]), even when the appropriate 
condition for the parametric resonance is satisfied. 
Besides, there is a coupling between the fast and 
slow magneto-sonic waves due to the open nature 
of the system, as is manifested from the terms on 
the RHS of Eqs. (I2S])-(I311). This coupling is very 
weak when the plasma /? is very large and the fre- 
quencies of the waves are very separated. Nonethe- 
less, this coupling leads to the excitation of slow 
(cf. panel Ci in Fig. [T|) and fast (cf. panel B2 in 
Fig. [T]) "noises" of very small amplitude. 

3. The situation is just opposite in the case of low 
plasma-/? (^^ ^ 1). In this case only the slow 
magneto-sonic wave frequency oscillates with the 
significant amplitude, while the variability of the 
fast magneto-sonic wave frequency is vanishingly 
small. Hence, only slow waves can be amplified (cf. 
panel C2 in Fig. ^ by the parametric action of the 
forced entropy oscillations and the fast magneto- 
sonic modes do not feel (in the sense of parametric 
action) these oscillations effectively (cf. panel Bi 
in Fig. ^ with the constant Alfven speed consid- 
ered in current work. In addition, we again observe 
the excitation of slow magneto-sonic (see panel Ci 
in Fig. ^ and fast magneto-sonic (see panel B2 in 
Fig. [2) "noises" due to the weak mode coupling. 

4. When (^^ ~ 1), we can distinguish two different 
situations: 
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(i) In the case when both fast and slow magneto- 
sonic waves propagate almost parallel to the mag- 
netic field Kx <C 1, the characteristic frequencies 
of the different modes are very close to each other. 
Therefore, there is a strong wave coupling. As a 
result, there is no possibility for only one kind of 
wave mode to exist, i.e. the presence of one of the 
wave modes immediately results in the excitation 
of the other. In this case, both kinds of waves can 
be involved together in the parametric resonance as 
shown in Fig. [3l 

(ii) However, when the waves propagate in the 
oblique direction with respect to the applied mag- 
netic field, the frequencies become somewhat more 
separated and then the coupling effect weakens sig- 
nificantly (Fig. m. 

The results obtained in this work can be applied un- 
der terrestrial laboratory conditions such as in acoustic 
investigations but the model also appears useful for the 
description of some astrophysical processes that are pe- 
riodic in time, as in objects with a variable intensity of 
entropy production and losses (as for instance in vari- 
able stars). However, at the present stage we aim only 
to demonstrate the basic physical properties of the para- 
metric amplification of waves and their mutual coupling 
in nonequilibrium plasmas. Despite of the apparent sim- 
plifications our new results are believed to be of general 
importance and they may have far reaching consequences 
for many applications. Here follow some possibilities: 

• Previous investigations of processes based on para- 
metric resonance in neutral and in conducting fluids 
have been carried out. These studies have revealed 
the possibility of an effective parametric amplifica- 
tion of sound waves either due to external mechan- 
ical periodic forcing or to mutual nonlinear interac- 
tions between the different harmonics of the waves 
themselves. Along with the references on this mat- 
ter given in the previous sections, other examples 
include [H, [l^, [2^. The model outlined in the 
present work is a unified approach to parametric 
amplification of sound and of compressional MHD 
waves in considered two dimensional systems. The 
generalization to 3D would also enable an experi- 
mental treatment of the parametric resonance trig- 
gered by an externally driven, oscillating entropy 
(or temperature). Thus, it would yield a valuable 
contribution to the rich spectrum of theoretical and 
experimental studies (for instance, see [IB, [1^ and 
many others) in the field of parametric instabilities 
in fusion plasmas. 

• The spectrum of astrophysical processes, where the 
considered mechanism of wave amplification and of 
wave coupling may be useful is even more rich. We 
have already given (see above) several references 
concerning the swing (parametric) interactions be- 
tween different MHD waves in equilibrium mag- 
netic structures, which directly have been applied 



to the magnetized solar atmosphere. On the other 
hand, there are many studies, which reveal the im- 
portance of parametric resonances at galactic and 
cosmological scales (for instance, see [5l[8l,[T8l[37|). 
It should be emphasized that the results obtained 
in this work may explain the properties of oscil- 
latory phenomena in astrophysical objects being 
away from thermodynamic equilibrium and sus- 
taining forced entropy oscillations. The entire va- 
riety of such thermally variable objects are usually 
described within the framework of thermal relax- 
ation oscillation theory [s^]- We give here refer- 
ences to types of the variable stars which mani- 
fest the mentioned oscillatory behavior: Thermal 
oscillation during the carbon burning and He- 
lium burning [10] in the stellar core; Also thermal 
relaxation in the horizontal-branch stars [l7j and 
those with convective cores (31|. In addition, our 
approach appears interesting for models of globally 
oscillating coronas [l^, although such a consider- 
ation would require a generalization of the current 
model for the case when the variability of the den- 
sity is also taken into account. In such systems the 
parametric amplification of the physical quantities 
could even become observable if it would be pos- 
sible to observe two, mutually related oscillatory 
processes with the frequencies related as toi « 20^2. 

The model presented here, of course, requires some 
further development, e.g. to enable the study of similar 
effects when the background density and velocity are also 
variable in time. In addition, the terms corresponding to 
a perturbation of the source term in the entropy balance 
equation, due to the propagation of waves, could also be 
taken into account (in the current model we neglected 
this term). In that case, additional kinds of modes can 
arise in the system such as thermal waves and aperiodic 
vortices which would enrich the spectrum of possible in- 
teractions. 
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APPENDIX: RELATION OF THE HEAT 
CONDUCTION WITH THE ENTROPY 
PRODUCTION AND ITS CURRENT 

The current derivation is based on standard linear re- 
sponse theory for weakly non-equilibrium systems e.g. see 
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01 • Consider the source term A AT in Eq. ((T2)) . It is well 
known that the coefficient of thermal conduction A arises 
here due to the assumption that it depends only on the 
overall equilibrium temperature A — Lgq/T^Q, where Lqq 
is a constant parameter. Setting Lqq = in the con- 
sidered case reduces the system to the equilibrium. On 
the other hand, the assumption that the coefficient A is a 
constant is known as the Fourier approximation, leading 
to the Fourier law of heat transfer (fT2|) . From a gen- 
eral point of view the equation which governs the heat 
transport reads as: 

- (^^^°) = (^-^^) - -^-J- 

which says that the changes in the local entropy are 
caused by the heat transfer through the system (if there 
is no source of heat within it). The divergence of the heat 
current on the r.h.s. can be reorganized as follows: 

-V-j, = ro(-V-j,+a), (A.l) 

which recovers the expression ([T]) with the microscopic 
entropy current js = jg/Tp and local entropy production 

<T=j,-v(^^). (A.2) 

These last equalities reveal the relation between the coef- 
ficient of thermal conduction A, and the entropy current 
and its production. As a matter of fact, the presence 
of the coefficient A in Eq. manifests the presence of 
both entropy production and its transfer, thus determin- 
ing the system to be away form thermodynamic equilib- 
rium. 



Further, in the Fourier approximation one can write: 
ao « ^ (VTo)' (A.3) 



-V.j,«— — (VTof (A.4) 

These expression can be written explicitly for the exam- 
ple dUl): 

do ~ -T-^f exp {2 fx) exp(2iwt), (A.5) 

Lqq 



-V • « — exp {2fx) exp(2«cji)+ 



^ ^pTm exp {fx) exp(wt) 
Tqo + Toi exp {fx) exp(iwi) ' 

where. 

As it has been mentioned in the text we omit the co- 
ordinate dependence of these quantities implying that 
exp (/a;) « exp(/xo), where Xq is a constant. This last 
assumption immediately leads to an equation of type (jlOp 
with a constant amplitude ai, which is under considera- 
tion here. 
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